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Abstract 

We discuss the possibility that gamma-ray bursts result from in- 
ternal shocks in an ultra-relativistic matter. Using a simple model we 
calculate the temporal structure and we estimate the efficiency of this 
process. In this model the ultra-relativistic matter flow is represented 
by a succession of shells with random values of the Lorentz factor. We 
calculate the shocks that take place between those shells and we esti- 
mate the resulting emission. Internal shocks can produce the highly 
variable temporal structure observed in most of the bursts provided 
that the source emitting the relativistic flow is highly variable. The 
observed peaks are in almost one to one correlation to the activity of 
the emitting source. A large fraction of the kinetic energy is converted 
to radiation. The most efficient case is when an inner engine produces 
shells with comparable energy but very different Lorentz factors. It 
also gives the most preferable temporal structure. 

1 introduction 

The Burst and Transient Source Experiment (BATSE) on Compton-GRO 
has revolutionized our ideas about the location of gamma-ray bursts(GRBs). 
The isotropy of the events and the paucity of week bursts strongly support 
the cosmological origin of GRBs (Meegan et al. 1992; Paczynski 1992; Piran 
1992; Nemiroff et al. 1994). The most fascinating mechanisms for produc- 
ing cosmological GRBs is decelerations of ultra relativistic matter having a 
Lorentz factor 7 > 100. This provides the only known solution to the com- 
pactness problem (Piran 1996). The kinetic energy of the ultra-relativistic 
matter is converted into internal energy by relativistic shocks. These shocks 
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can be due to the inter stellar medium (ISM), "external shocks" or inside 
the shell itself due to non-uniform velocity "internal shocks" . Electrons are 
heated by the shocks and the internal energy is then radiated via synchrotron 
emission and inverse Compton scattering. Though the ultra relativistic mat- 
ter flow was originally considered in the dynamical context of fireball(Shemi 
& Piran 1990; Paczyhski 1990), one can imagine GRB models in which the 
fireball is replaced by some unknown non-thermal acceleration mechanism 
but the radiation still originate from slowing down of the ultra relativistic 
matter (Piran 1996). 

Most of bursts have a highly variable temporal profile with a time scale 
of variability 5T significantly shorter than the overall duration T. Typical 
value is ST/T ~ 10~^. Clearly any GRB model must be able to explain this 
complex temporal structure. 

The naive external shocks scenario requires only one uniform shell. The 
complex structure is due to surrounding matter (ISM in some models and 
star light in others (Shaviv & Dar 1995)). Two of us have shown that the 
efficiency of this process is less than 1% (Sari & Piran 1997a). Using the 
observed variability scale 6T one can estimate the typical size of a single 
emitter. The number of emitters must be less than the number of peaks. The 
total area of all emitters is therefore only a ~ 1% of the shell's area. Even 
a modified version of this model, in which the shell is composed of many 
smaller shells that collide one after the other with the ISM is impossible. 
Most of the collisions would occur at a radius which satisfies R > A'j'^ 
where A is the total width of the shellsQ At this radius angular spreading 
(Katz 1994; Fenimore et al. 1996) will cause each peak to be spread on 
time scale of R/'j'^c = A/c which is also the total duration of the bursts. 
The observed bursts will therefor be smooth. The only way out, in the 
context of the external shock model is if the shell is replaced by a narrow 
and highly variable beam. However, the angular width of this beam bust be 
than 5T/iT-f) ^ 10"^! 

The obvious alternative is the internal shock scenario (Rees & Meszaros, 
1994; Narayan, Paczyhski & Piran, 1992). These shocks could occur within 
a variable relativistic wind produced by a highly variable source. Internal 

^Distances, time, velocities /3 and and the corresponding Lorentz factors 7, are mea- 
sured in the observer's rest frame. Thermodynamic quantities are measured in the fluid 
rest frame. 
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shocks do not suffer from the angular spreading problem. A fast shell that 
was injected after a slower one will eventually catch up and collide with it. 
If both shells have Lorentz factor of order 7, the time until collisions will be 
of order of 7^L/c, where L is the initial separation between the shells. The 
collision will therefor occur at radius R = '-f'^L and the angular spreading 
time is of order of L/c. Since the total duration is given by the total width 
of the shells A/c ^ L/c, the observed peak width is considerably smaller 
than the total width and complex temporal structure is possible. In view of 
the difficulties within the external shock scenario it is worthwhile to consider 
the question whether internal shocks could actually produce the temporal 
structure observed in GRBs and what would be the efficiency of this process. 

To examine the energy conversion in internal shock we construct, fol- 
lowing (Mochkovitch Maitia & Marques, 1995) a simple model in which the 
ultra-relativistic wind is represented by many shells having a random distri- 
bution of intrinsic parameters. We calculate the energy conversion due to 
shocks that form between the shells, considering at each time a binary (two 
shells) encounter. We show that internal shocks can produce the observed 
highly variable temporal profiles. We reproduce the result of Mochkovitch, 
Maitia & Marques (1995) who concluded that the efficiency of this process 
is low (< 10%) if the spread in Lorentz factors is small. However, we show 
that higher efficiency could be achieved if the spread in the Lorentz factors 
of the ultra-relativistic matter is larger. 

We explore in section 2 the basic unit of our model, the interaction of two 
shells. We explain in section 3 the algorithm for evolution of multiple shells. 
In section 4 and 5, we discuss the temporal structure and the efficiency, 
respectively. Finally we discuss the implication to observation in section 6. 

2 Two shell interaction 

Internal shocks arise in a relativistic wind with a non uniform Lorentz factor 
and convert a portion of the kinetic energy to a radiation. We represent the 
irregular wind by a succession of relativistic shells. A collision of two shells 
is the elementary process in our model. 

A rapid shell (denoted by the subscript r) catches up a slower one (s) and 
the two merge to form a single one (m) . The system behaves like an inelastic 
collision between two masses and m^. Using conservation of energy and 
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momentum we calculate the Lorentz factor of the merged shell to be 



(1) 



where 7i(^ 1) and mi are the Lorentz factors and the masses of these shells. 
The internal energy of the merged shell is the difference of kinetic energy 
before and after the collision: 



Eint = mrC^{-fr - 7m) + msC^{% " 7m)- 



(2) 



The efficiency of conversion of the kinetic energy into the internal energy in 
a single collision is given by 



1 - {nir + ms)jm/{mrjr + "^s7s)- 



(3) 



The emitted radiation will be observed as a pulse with a width 6T. Three 
time scales, the cooling time, the hydrodynamic time, and the angular spread- 
ing time determine 6T Q. The internal energy is radiated via synchrotron 
emission and inverse Compton scattering. In most of cases, the cooling time 
scale of electrons is much shorter than the hydrodynamic time scale (Sari et 
al. 1996; Sari & Piran 1997b), so we can neglect the cooling time. 

The hydrodynamic time scale is the time that the shock crosses the shell. 
In fact there are two shocks as the interaction between the two shells takes 
place in the form of two shocks: a forward shock and a reverse shock. Using 
conservation of mass, energy and momentum at the shocks and equality of 
pressures and velocities along the contact discontinuity, one can derive the 
Lorentz factors of the forward and the reverse shocks 7/s,7rs (Sari & Piran 
1995): 
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^Clearly, if the source is located at a high redshift z, all the observed time scales are 
stretched by a factor I + z. However, as we are interested mainly in the ration between 
the width of the peaks and the overall duration and in the efficiency energy conversion 
this factor is unimportant. 
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For simplicity, we estimate the emission time scale by the time that the 
reverse shock crosses the rapid shell: 



Ste = lr/c{Pr - Prs)- (5) 

where is a width of the rapid shell. The emitting region moves towards 
the observer, the observed time scale is shorter by a factor of 1/27^. If both 
shells have Lorentz factor of order 7, this observed time scale is of order /^/c. 

If the coUision of the shells takes place at a large radius R, angular spread- 
ing effects the width of the pulse. The separation of the shells is L, the effect 
of angular spreading on the pulse width is ~ L/c as we mentioned in sec- 
tion 1. If the separation L is larger than the width of the shells /, the pulse 
width 6T is determined by angular spreading. The shape of the pulse be- 
come asymmetric with a fast rise and a slower decline (Fig. 1) which GRBs 
typically show. The peak amplitude is given by Eintc/l. If angular spreading 
is significant, the amplitude is lower by a factor of 1 — (1 + l/L)~'^. The 
observed luminosity, C is given by 

r 0, {t< 0), 

C{t) = {h[l- 1/(1 + 2^lctlRf] , (0 < t < 5te/27^), 

[ h [1/(1 + {2^lt - 5Qc/Rf - 1/(1 + 2^lct/Ry] , (t > 6tj2^l), 

(6) 

where h = Eint2^l,/ 5te. 

We assume that the internal energy produced by the collision is radiated 
isotropically in the shell's rest frame. We also assume that all the internal 
energy is radiated rapidly. In this case we obtain at the end of the collision 
a single cold shell whose Lorentz factor is 7^, given by Eq. |l|. The shocks 
compress the initial shells and it results in a thinner shell with a width Z^: 

The density is discontinuous through the contact discontinuity. For simplicity 
we average the density: 

Pr^r^Jr ~l~ Ps^s^s /n\ 
Pm = -J • (8) 

We describe the resulting shell as a homogeneous shell with a constant density 
Pm and a Lorentz factor 7^. 
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While the shells propagate, their density decrease due to the spherical 
geometry. However the density ratio between any two shells remain constant 
and the process depends just on the density ratio and not on the absolute 
value of the density. We expect therefor that this one dimensional model 
captures the basic features of realistic events. 

3 The Multiple shell model 

We represent the irregular wind by a succession of relativistic shells with a 
random distribution of Lorentz factors and densities. A collision between two 
shells produce shocks which convert some of the kinetic energy of the shells 
to thermal energy. The thermal energy is then radiated via synchrotron 
emission and inverse Compton scattering. In section 2, we calculated the 
observed radiation produced by a collision between two shells. For multiple 
shells, numerous coUisions take place. Each coUision produces a pulse similar 
to that shown in Fig. 1. To construct the temporal structure in a given model 
we calculate the time sequence of the two shell collisions and we superimpose 
the resulting pulses from each collision. 

Consider a wind consisting of N shells. We assign an index i, {i = 1, N) 
to each shell according to the order of the emission from the inner engine. 
Each shell is characterized by four variables: a Lorentz factor 7j, a density 
Pi, (the shells are cold and therefore the pressure satisfies Pi = 0.) a width 
li and the time, tj, when the shell was ejected from the source. We choose 
the origin of time such that the last shell - the N-th shell - is emitted at 
fjv = 0. Clearly all other shells are emitted earlier so the tj's are negative 
ior i < N. At t — the shells are at the following "initial" positions (inner 
edge): Ri = —iiPiC. We denote the distance at this stage between the shell 
i and the shell i + 1 by Lf Li = Ri — -Rj+i — /j+i. The time li/c corresponds 
to the period that the "inner engine" has operated, while Li/c corresponds 
to the period for which the "inner engine" was quiet. The size of the inner 
engine should be, of course, smaller than / to produce a shell with a width 

We turn now to evaluate the evolution of this multiple shell system. There 
will be numerous collisions between different shells. We denote by the index 
j the j-th collision which takes place at a time tj and position Rcitj)- For 
simplicity we define a zeroth collision at the initial time to — Q this enables 
us to put the first collision on par with the rest. Given the velocities of the 
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shells f3i and their separation Li at tj_i we calculate for all pairs + 1) 
satisfying /Jj+i > /3j the collision time: 

6ti^i+i = Li/c{pi+i - Pi). (9) 

We then find the minimal collision time among all the 5tiS: 

Stj = mm[dti^i+i]. (10) 

Let s and s + 1 be the index of the shells for which this minimum takes place. 
The j-th collision is between the s and s + 1 shells and takes place at the 
time 

tj = tj.i + 5tj (11) 

and at the position 

Rc(tj) = Rs+i(tj) = Rs+i(tj^i) + c(3s+iStj. (12) 

An observer at Ro away from the central source will begin to detect radiation 
from this collision at a time 

tobs,s={Ro-Rc{tj))/c + tj. (13) 

Note that, for reasons that will be clear later we give this observed time the 
index s of the outer shell that participated in the collision. At the end of 
the computation we set the minimum of tobs,s to be the origin of the observer 
time as this is when the first radiation is observed. We assign the pulse shape 
presented in section 2 beginning at this time to the overall observed signal. 

We rearrange the shells after each coUision as follows. Each shell moves 
from its earlier position Ri{tj-i) to 

Ri{tj) = Riitj.i) + cpiStj (14) 

For the first collision tj_i = and we use the "initial positions" here. All 
the shells except the s shell and the s + 1 shell keep their Lorentz factor, 
the density and the width: 7^, pi and li. The s and the s + 1 shell merge to 
form a single shell which we label s + 1. We delete the s shell and we will 
skip from now on this index when it will be encountered. Using the results 
of section 2 we calculate now the new parameters for the s + 1 shell: js+i, 
ps+i and Is+i- 
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We return now to the calculation of the distances between shells: Li and 
the corresponding collision time 5ti. We find the next collision and proceed 
until there are no more collisions, i.e., until all the shells have merged to 
form a single shell or until the shells are ordered with increasing values of 
the Lorentz factor. 

In order to rcdiicc the number of parameters describing the wind we adopt 
the following simplifications. We assume a constant initial width /j = / and 
constant initial separation Lj = L. In other words, we assume that the "inner 
engine" operates for a fixed period I / c and then it is quiet for a fixed period 
L/c. 

We assume that the Lorentz factor of each shell is uniformly distributed 
between 7^^^ and 7^03; • For the distribution of the density we have two kinds 
of models. In the first kind we consider models in which the density, the mass 
or the energy of each shell is a random variable which is uniformly distributed 
between 1 and some maximal value X^ax- In the second kind the density 
is correlated with the Lorentz factor as pi oc 7,f~^ so that the mass and 
kinetic energy of each shell are proportional to 7J' and 7^^^ respectively.. For 
T] — 1 the shells have equal density, for 77 = the shells have equal mass and 
for ?7 = — 1 the shells have equal energy. We will see later that the random 
density/mass/encrgy models give almost same efficiency as the corresponding 
constant model, and the dependence on Xmax is low. We therefore present 
most of the result for the correlation models which are more simple. 

Any time scale in this system is proportional to a sum of L and Z, the 
amplitudes of observed peaks depends neither on the absolute value of I nor 
that of L. If both of L and / are transformed by the same factor, we will 
get similar temporal profile. We are interested in neither the absolute value 
of the duration nor that of the peak, then we can use a single parameter 
L/l instead of the two parameters L and I. The efficiency of the energy 
conversion by the shocks is independent ol L/l. The total energy emitted is 
the sum of the energy emitted in the elementary two shell collisions which 
are independent of L/l. The density p is the only quantity in this system 
with dimensions of mass. Therefore only density ratio rather than density 
can determine the time scale, and the overall normalization of the density 
p is also unimportant. Consequently we find that there are only five initial 
parameters: A^, L/l, 'Jmax, Imin and either rj (if the density is correlated with 
the Lorentz factor) or X^nax (for random distribution of p,m or E). 
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4 Temporal Structure 



The calculated temporal structure is a superposition of the pulses from the el- 
ementary two shell collisions. Several typical temporal profiles are presented 
in Figs. 2a-2f. Fig. 2a depicts the luminosity for a model with = 100 and 
L// = 5 as a function of the observed time for shells with a constant energy 
(?7 = —1). During the evolution of the N shells, almost N collisions happen 
and almost N pulses are produced. Though some of the peaks are not ob- 
served because of their higher neighbors. The number of peaks is given by 
~ A^. In Fig. 2b the luminosity is plotted with the same parameters as Fig. 
2a except that N = 20. Fig 2a and 2b show that practically determines 
the number of peaks. The width of a given peak ST is determined by the 
emission time scale ~ Z/c or by the angular spreading L/c. If the value 
of the separation L is larger than that of the width ST/T is ~ 1/N and 
independent of L/l. In Fig. 2c, the luminosity is plotted with the same pa- 
rameters as Fig. 2b except that L/l = 1. The profile is almost independent 



The different amplitude of the peaks originate mainly from the difference 

in Eint, since all peaks have width of the same order of magnitude L/c. 
For given random Lorentz factors, Eint depends on the index 7] for models 
with density correlated to the Lorentz factor. The random Lorentz factors 
are uniformly distributed between 7^,^ and -^raax- If Imax/lmin > 1, the 
ratio 7r/7s is typically large. In such a case, we obtain a simple formula. 



For ?7 = — 1, Eint takes almost the same value for most of the collisions so 
most of the peaks have comparable amplitude (see Fig. 2a). Thus, for rj = —1 
we practically observe the "shell structure" produced by the inner source as 
almost all two shell collisions produce an observable peak. If on the other 
hand rj — 1 (see Fig. 2d) a few collisions produce significantly higher peaks 
than others and only those are observed. Thus the number of the observed 
peaks is much less than N, the observed temporal structure corresponds only 
weakly to the activity of the inner source. Visual inspection suggests that 
the temporal structure produced by 77 = —1 resembles better the observed 
temporal structure. 



of L/l. 




iv = i) 
iv = - 



(15) 
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Fig. 2e and Fig. 2f present the luminosity profiles for model with random 
energy and random density respectively.. Both use Xmax = 1000 and both 
with spread of Lorentz factor ^max/lmin — 10- Again the random density 
is more spiky than the random energy case. However as could be guessed 
intuitively, the random energy and the random density cases are more spiky 
than the constant energy (77 = —1, Fig. 2a) and constant density (77 = 1, 
Fig. 2d) respectively.. 

The total width of the wind is A = A^Z + (A^ - l)L. Neighboring shells 
collide on a time scale of -^7^/ c, with fluctuations of the same order of magni- 
tude due to the fluctuations of 7. The matter is moving toward the observer, 
the resulting observed time scale is L/c. On the other hand, the difference 
in observed time due to the location of the given shell within the wind is a 
order of A/c. Since by definition A ^ L we observe pulses arising from col- 
lision between shells mostly according to their positions inside the wind, i.e., 
according to the time that those shells were emitted by the inner engine. In 
Fig. 3, we plot the time that radiation is observed from a shell tobs.s versus 
the time that the shell was ejected from the source ig- We associated the 
emission from a given collision with the slower shell. The clear correlation 
between the observed time and the time that the source ejected the specific 
shell shows that the observed temporal structure reflects the activity of the 
source. This conclusion is valid even if there is a large spread in 7. 



5 Efficiency 

We have shown that internal shocks could produce the highly variable tem- 
poral structure. We ask now what is the efficiency of this process? The 
relativistic shells collide with each other and merge into more massive shells. 
The overall efficiency of conversion of kinetic energy in internal shocks can 
be calculated from the initial and final kinetic energies as: 

where the superscript (/) and (i) represent the initial and final values, re- 
spectively. 

This efficiency depends on the parameters of the model A^, '~fmin, Imax 
and 1] (or Xmax) and on the specific realization: the set of random Lorentz 
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factors that are assigned to each shell. For each choice of the parameters of 
the model, we have evaluates the efficiency for 100 realizations. The mean 
efficiency and its standard deviation are listed in Table 1. 

The efficiency is only a few percents if the spread in the Lorentz factor 
is relatively low (a factor of 2-3). This agrees with results of Mochkovitch, 
Maitia and Marques (1995) who concluded that the efficiency of this process 
is less than 10%. However, higher efficiency could be reached if the spread 
in 7 is larger. The spread required, for example, for 10% is ^max/lmin ~ 6 
for N = 100 and r] = -1. 

The efficiency is independent of L/l, which is important only for the 
temporal structure. For models with density correlated with the Lorentz 
factor the relevant parameters that determine the efficiency are N, rj, •jrnin 
and jmax (see Table la). For the random density models, we evaluated the 
efficiency with X^ax — 1 to 1000. The efficiency is almost independent of 
Xmax- It can be seen that the efficiencies of the models with random n, m 
and E are very similar to those with 77 = 1,0 and —1, respectively (Table 



As the number of shells, A^, increases, the kinetic energy of relative mo- 
tion of the final mergers becomes negligible. The efficiency approaches an 
asymptotic value. Further if we consider large values of 7mm, the efficiency 
depends only on the ratio ^max/lmin reaching an asymptotic value for large 
Imax/lmin (Fig- 4). For models with correlation between the density and 
the Lorentz factor, this asymptotic value depends only on the index 77. In 
the limits — > ±00, a single shell carries all the energy of the system and 
the energy and momenta of the other shells are negligible. Therefore, the 
massive shell is unaffected by collisions with the lighter ones. We expect a 
low efficiency in this limit (Fig. 5). 

In order to understand the dependence of the efficiency on rj we consider 
a simphfied model. Let all the shell coUide and merge into a single shell and 
only then emit the thermal energy as the radiation. Using conservation of 
energy and momentum we can calculate the Lorentz factor and the efficiency: 



Averaging over the random variables 7^, and assuming a large number of 



lb). 




(17) 




(18) 
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shells N ^ oo we obtain: 



fry . )V+i _ 1 

/ \ -I V I max I mm J J- 

(e) ~ 1 



1] ~\~ ^ ^ (^(Sfmax / 'ymin)''^ 1) ((Tmaa'/O'mm)'^''' 1) 

(19) 

This formula explains qualitatively the behavior of our numerical results. A 
comparison between the simulations in which the radiations is emitted right 
after each collision and the simplified single merge model is presented in Fig. 
5b. 

An interesting point to notice is that for many shell collision (both in the 
numerical simulation and the simple model Eq. ^) the most efficient case is 
where the shells have the same energy while for the two shell collision (Eq. |^) 
the highest efficiency is for shells with equal mass rather than equal energy. 



6 Conclusions 

Using a simple model, we have shown that internal shocks can produce the 
highly variable profile observed in most GRBs. There is a strong correla- 
tion between the time at which we observe a pulse and the emission time 
of the corresponding shell from the inner engine. This correlation persists 
even when there is a large spread in the Lorentz factor. Thus the observed 
temporal structure reproduces the activity of the source. 

We have shown that the number of peaks is almost the same as the 
number of shells that the inner engine emitted. The separation between the 
peaks corresponds to the duration for which the inner engine was quiet. The 
variability of peaks height can tell as wether all shells have comparable energy 
(low variability in observed peak heights) or not. A systematic statistical 
study comparing the temporal structure produced by internal shocks with 
observations will be carried out in a subsequent paper. 

The efficiency of this process is low (less than 2%) if the initial spread in 

7 is only a factor of two. However the efficiency could be much higher. The 
most efficient case is when the inner engine produces shells with comparable 
energy but with very different Lorentz factors. In this case [i] = —1, and 
spread of Lorentz factor 'jmax/lmin > 10^) the efficiency is as high as 40%. 
For a moderate spread of Lorentz factor jmax/lmin = 10, with rj = —1, the 
efficiency if 20%. 
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Tabid a 



N 






Ijnax 


pflfiriPTirv [ynl 


100 


1 


20 


1000 


10.9 ± 1.4 


100 


1 


50 


1000 


10.0 ± 1.3 


100 


1 


100 


1000 


8.5 ± 1.1 


100 


1 


200 


1000 


6.1 ± 0.7 


100 


1 


500 


1000 


1.7 ±0.2 


100 





20 


1000 


19.2 ± 2.9 


100 





50 


1000 


16.0 ±2.3 


100 





100 


1000 


12.4 ± 1.6 


100 





200 


1000 


7.7 ±0.9 


100 





500 


1000 


1.8 ±0.2 


100 


-1 


20 


1000 


25.0 ±3.5 


100 


-1 


50 


1000 


19.5 ±2.4 


100 


-1 


100 


1000 


14.1 ±1.6 


100 


-1 


200 


1000 


8.2 ±0.9 


100 


-1 


•lUO 


1000 


1.8 ± 0.2 


20 


1 


100 


1000 


7.6 ±2.5 


20 


-1 


100 


1000 


11.0 ±3.3 


1000 


1 


100 


1000 


8.7 ±0.3 


1000 


-1 


100 


1000 


15.1 ±0.5 


1000 


-1 


10 


10000 


39.3 ±2.6 



Table. 1 b 



iV 






'y 771 ax 




efficiency [%] 


100 


random n 


100 


1000 


1000 


8.2 ± 1.2 


100 


random m 


100 


1000 


1000 


12.2 ± 1.6 


100 


random E 


100 


1000 


1000 


14.2 ± 1.6 



Mean efficiency and standard deviation of the efficiency for 100 realizations of different 

models. 
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Figure captions 



Fig. 1 A peak produced by a collision between two shell. The luminosity 
plotted versus the observer's time. The solid line corresponds to R = c6te 
and the dotted line corresponds to i? = 0. 

Fig. 2 The observed temporal structures. The luminosity versus the 
observer's time, for different models 

a: -fmin = 100, -f^ax = 1000, N = 100, r] = -l and L/l = 5 
b: 7™„ = 100, -/max = 1000, iV = 20 , r/ = -1 and L// = 5 
c: = 100, -f^ax = 1000, N = 20 ,ri = -1 and L/l = 1 
d: 7™„ = 100, -fmax = 1000, N = 100, r/ = 1 and L/l = 5 
e: -fmin = 100, -fmax = 1000, = 100, 

random energy with Emax = 1000 and L/l = 5 

f: -fmin = 100, Jmax = 1000, = 100, 

random density with Pmax = 1000 and L/l = 5 

Fig. 3 The time of ejection of a shell by the inner engine, tj versus the 
observed time of the photon produced in that shell, tobsj, for A^ = 100, 
7mm = 10, jmax = 1000, T] = —1 and L/l = 5. The initial position and 100 
correspond to the inner and outer edge of the wind. 

Fig. 4 The efficiency versus 'ymax/lmin with la error bars of 100 random 
simulations, for ■ymin = 10, N = 500 and r] = —1. 

Fig. 5 The efficiency plotted versus the index r], for 'jrnin = 100, A^ = 100. 
The solid line corresponds to 'jmax/lmin = 1000, the dotted line corresponds 
to ^max/lmin = 100 and the dashed line corresponds to 'jmax/'lmin = 10. a: 
Numerical simulation, b: Approximation by analytic formula (|19D 
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